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The thermal conductivity k of the Spin- Peierls (SP) com- 
pound CuGeOs was measured in magnetic fields up to 16 T. 
Above the SP transition, the heat transport due to spin ex- 
citations causes a peak at ^22 K, while below the transition 
the spin excitations rapidly diminish and the heat transport is 
dominated by phonons; however, the main scattering process 
of the phonons is with spin excitations, which demonstrates 
itself in an unusual peak in k at ^5.5 K. This low-temperature 
peak is strongly suppressed with magnetic fields in excess of 
12.5 T. 

PACS numbers: 66.70.+f, 75.30.Kz, 75.50.Ee 



The recent discovery of the inorganic Spin-Peierls 
(SP) compound CuGeOs Q has revived interests in 
this phenomenon Spin-Peierls transition occurs in 

quasi-one-dimensional magnetic systems when the energy 
gained by splitting the degeneracy of the original mag- 
netic ground state exceeds the lattice deformation energy. 
The SP ground state is therefore a nonmagnetic singlet 
state with a magnetic energy gap between the singlet and 
the lowest excited triplet states. The appearance of the 
magnetic energy gap and the accompanying dimerization 
of the 5=1/2 Cu 2+ ions have been confirmed by magnetic 
susceptibilities |Q, inelastic and elastic neutron scatter- 
ings X-ray diffractions ||, and other experimental 
probes. The size of the spin gap has been found to be 
A ~23 K ||. This gap can be suppressed with mag- 
netic fields; above a critical field H c c^l2.5 T, the system 
undergoes a first-order transition to an incommensurate 
(I) phase (|]. The incommensurate lattice pattern has 
been measured by X-ray experiments and interpreted as 
soliton lattice structure 

Although CuGeOs shows clear evidence for SP transi- 
tion, there has been accumulating evidence that CuGeOs 
is not a prototype SP system. For example, inelastic- 
neutron-scattering [[| find the nearest-neighbor (NN) in- 
teraction J ~120 K and a relatively strong interchain 
coupling (Jft ~ 0.1 J). Static susceptibility, which should 
peak at k^T ~ J for an ordinary SP system, shows a 
broad maximum centered at a temperature much smaller 
than J. To describe the magnetic properties of this sys- 
tem, a model which includes next-to-nearest neighbor 
(NNN) interaction a J has been proposed [ p!o|JTT| . Since 
both NN and NNN interactions have the same sign, large 
a is a source for frustration. For a < a c (a c estimated 
as 0.241 [|ll|) the ground state of the magnetic system 



remains gapless, while for a > a c the spectrum of the 
spin waves becomes gapped even above the SP transition 
temperature T$p. Recent estimates of the frustration 
parameter suggests a to be as high as 0.36 JlC|] . 

While the phase diagram of CuGeOs is now well estab- 
lished, no transport measurements have been reported to 
date on this system. Such measurements are important 
for the determination of the basic scattering mechanisms 
in this system, both above and below the SP transition 
and in the high-field I phase. Since the SP transition 
is basically a magnetoelastic transition, heat transport 
which is expected to take place by phonons and spin ex- 
citations is ideal to probe these quasiparticles. 

In this Letter we report thermal conductivity mea- 
surements of high quality CuGeOs single crystals along 
the chain direction Q in magnetic fields up to 16 T. 
One of our main results is an observation of the spin- 
excitation continuum above the SP transition causing a 
broad maximum at ~22 K in the thermal conductivity 
k(T). This spinon heat transport disappears rapidly be- 
low Ts-p, causing a kink in ft(T), and phonons become 
the dominant heat carrier. The most intriguing observa- 
tion is that the phonon thermal conductivity below Tsp 
shows an unusual peak at ~5.5 K. This peak is strongly 
suppressed with magnetic field above H c ^12.5 T where 
the SP order disappears. In the I phase above H c the 
thermal conductivity increases, which supports the idea 
of a new spectrum of magnetic excitations that are mo- 
bile and can carry heat through the crystal. 

The thermal conductivity results are supplemented 
with specific heat and magnetic susceptibility measure- 
ments all done on the same samples. The specific heat 
data are useful for deriving the thermal diffusivity (and 
thus the information on scattering times) from ft; in the 
analysis of the low-temperature specific heat data, we 
suggest a new fit that is consistent with the simple the- 
oretical idea of an almost mean-field SP transition and 
try to resolve some controversy over the analysis of the 
specific heat data [T^| . 

Single crystal fibers of CuGeOs were grown using the 
laser- heated pedestal growth technique, a miniaturized 
float zone process |1J]. Powders of pure CuGeOs are 
cold-pressed and sintered into a pellet, from which we 
cut several source rods of about 2.6 mm in diameter. 
Previously grown single crystal fibers are used as seeds. 
Our samples were grown using a 20 mm/h pull rate (seed) 
and about 5 mm/h feed rate (source rod) in a flowing O2 
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atmosphere. The as-grown crystals are ~4 cm long and 
1-2 mm in diameter. 

Thermal conductivity measurements were performed 
using a "two thermometer, one heater" method. Sam- 
ples of size 7 x 2 x 0.2 mm 3 were cut from the as-grown 
crystal, where the long dimension is along the chains (c- 
axis). The base of the sample was anchored to a copper 
block held at the desired temperature. A strain-gauge 
heater was used to heat the sample. A matched pair of 
microchip Cernox thermometers were carefully calibrated 
and then mounted on the sample. The maximum temper- 
ature difference between the sample and the block was 0.2 
K. Typical temperature gradient used between the two 
thermometers was ~0.5 K above 20 K and 0.2 K below 20 
K. Radiation shield was used and was kept at the average 
temperature of the sample. Specific heats were measured 
in the temperature range 1.7 - 20 K using the relaxation 
method. The addenda heat capacity was measured inde- 
pendently, with no sample attached, in separate runs. 
The particulars of our apparatus, which allows us an 
extremely high sensitivity measurements, are described 
in greater detail elsewhere |l5[ |. Magnetic susceptibility 
measurements were performed using a Quantum-Design 
SQUID magnetometer in the field range of - 7 T. 

Figure 1 shows the thermal conductivity of a CuGeOs 
single crystal in the temperature range 1.5 - 40 K. We 
also plotted in Fig.l the specific heat data measured on 
the same crystal in the temperature range 1.8 - 16 K. 
Note that the transition is clearly seen in the specific 
heat and its sharpness is similar to the best crystals re- 
ported to date. Before we turn to the discussion on the 
thermal conductivity, it is useful to analyze the specific 
heat data and see its correspondence to previously pub- 
lished results. Ignoring the transition region, we fit the 
data above the SP transition (Tsp= 14.08 K) to the form: 
c/T = 7 + /3T 2 , while at low temperatures we use the 
form: c/T = 5(A/k B T) n e- A / k * T + pT 2 . Note that we 
introduced the ratio {A/k B T) n as a prefactor to the ex- 
ponential decay of the spin contribution at low temper- 
atures. While all previous publications ]l6],|l7| used n= 
1, such a fit is erroneous and will hold only at very low 
temperatures, typically below 1% of the transition tem- 
perature. At higher temperatures an exact expression 
for c/T of the spin system is needed for a meaningful fit. 
In fact, we believe that the discrepancies in the litera- 
ture [|l3| between the value of (3 above and below Tsp 
stems from the use of incorrect fitting formula for the 
low temperature portion. Since the actual transition is 
3-dimensional in nature, the most reasonable model to 
use is that of BCS which gives n=5/2 (l8| . While the 
fit above the transition is straightforward, at low tem- 
peratures using the full formula with the gap function is 
problematic, for the fitting program will encounter a sin- 
gular matrix. Since the effect of the SP transition on the 
Debye temperature is expected to be negligible, a physi- 
cal approach is to use the value obtained for /3 above Ts p 
for the low temperature fit. 

Following the above discussion, we first fit the tempera- 



ture range above the transition (14.5 - 20 K), which gives 
7 = 105.5 mJ/mole-K 2 and = 0.2164 mJ/mole-K 4 . We 
then fit the low temperature portion (1.8 - 5 K) using the 
same (3 to obtain the gap. Using n=5/2 results in 5 = 
377.86 mJ/mole-K 2 and A = 25.7 K. This gives the ratio 
2A//cbX' c of 3.65, indicating a slightly stronger coupling 
than the BCS weak coupling limit of 3.52. It is useful 
to check the consistency of the parameter 5 by convert- 
ing it to the units of number of states per Cu atom, m. 
This is easily done by multiplying S by J (J as extracted 
from the literature lies in the range 80 to 150 K. Here we 
choose an averaged value of J ~115 K) and by the volume 
of a unit cell (a x b x c=1.2xl0" 22 cm 3 ). For the BCS 
formula with n=5/2 we get m=2.07 (note the \ / ^2tt factor 
in the prefactor of BCS formula p8[). This analysis sug- 
gests that the density of magnetic excitations is exactly 
2 states per Cu-atom at low temperatures. Thus, the use 
of the exact expression of c/T for the fit at T < Tsp 
leads to a totally consistent analysis with a reasonable 
value of A. 

We turn now to the thermal conductivity k in zero field 
shown in Fig.l. As the temperature is lowered from 40 K, 
k(T) increases and peaks at ~22 K below which it starts 
to drop rapidly. At the SP transition, k(T) shows a kink 
towards a faster drop. However, instead of continuously 
dropping (presumably exponentially), k(T) shows a min- 
imum and a subsequent increase to a new maximum at 
~5.5 K followed by a sharp decrease as the temperature 
is lowered further. As a first approach to understand 
the thermal conductivity we use a simple kinetic approx- 
imation. Assuming two types of heat carriers, phonons 
and spin excitations, the thermal conductivity can be 
written as: k = k s + k v = c s D s + c p D p . Here the sub- 
script s denotes the magnetic (spinon) component and 
p the phonon contribution, c is the specific heat and 
D the diffusivity of the relevant excitations. While the 
magnetic contribution to the specific heat could be fit- 
ted with a 7T term near T^p, the full description of the 
magnetic specific heat is more complicated and in fact 
is a matter of current controversy [ pj]Jl9| ]. For the pure 
one-dimensional system with only NN interaction, the 
specific heat will rise from low temperatures to a broad 
maximum at ~ J/2. This maximum shifts towards lower 
temperatures as the NNN interaction increases from zero. 
In fact, current models that fit the magnetic susceptibil- 
ity and Raman scattering data suggest a to be in the 
range 0.24 - 0.36 For example, Castilla et 

al. (IT) argued that for CuGeOs the competing exchange 
between NN and NNN interactions is moderately large 
but smaller than a c . Other authors report 0,|l^] that 
a larger value of a (~0.35) gives good fits to the experi- 
mental data. 

Since the magnetic contribution to the thermal conduc- 
tivity will be proportional to the magnetic specific heat, 
it is reasonable to assume that the peak at ~22 K is just 
a reflection of this magnetic density of states. The fact 
that this peak appears at a lower temperature than the 
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susceptibility (typically —55 K we also measured the 
susceptibility on the same sample confirming the peak at 
the same temperature) is a consequence of the tempera- 
ture dependence of the diffusivity D s — T~ q with q > 1. 
Assuming that k s is negligible by —100 K and approx- 
imating n p to be constant except at very low tempera- 
tures, we can subtract in the temperature range 15 to 50 
K a constant phonon background of k p c± 0.4 W/cmK. 
The remaining contribution approximates k s , which has 
a peak at —22 K. Shifting this peak to —50 K in c s re- 
quires a power q of 4 or 5 in the formula of the diffusivity. 
In fact, a simple scaling approach to estimate the spinon 
diffusivity above T SP yields D s - (J/h)a 2 (k B T/J)- q , 
where a is the distance between adjacent spins. Divid- 
ing the estimated k s by this expression for the diffusivity 
yields an approximate contribution of the magnetic ex- 
citations to the specific heat. Using <?=4, this magnetic 
specific heat shows a broad maximum at —50 K. Near 
Tsp the obtained curve is linear in T with a coefficient 
7* ~ 31 mJ/mole-K 2 (with J ~115 K), within a factor 
of 3 of the 7 estimated from the specific heat data near 
T$p. Since we omitted constants and numbers of order 
unity, it is indeed very encouraging that the two esti- 
mates are only a factor of 3 different from each other. 
Moreover, we subtracted a constant n p to obtain k 3 , a 
very simplified approximation. Thus, being able to ob- 
tain a reasonable value for the linear coefficient of the 
spinon specific heat give us confidence that indeed the 
maximum of the thermal conductivity at —22 K is of the 
same origin as the maximum in the susceptibilty. In ad- 
dition, the fact that this maximum is shifted so much 
to lower temperatures suggests that the mean free path 
of the magnetic excitations increases very rapidly as the 
temperature is lowered. 

The SP transition is manifested in n by a sharp down- 
ward change in slope. This drop is clearly related to the 
emergence of a spin gap and thus has the same origin as 
the one observed in the susceptibility (l) or the strong 
decrease of \jT\ below Tsp as seen in NMR measure- 
ments [pl[ . It is however difficult to determine the exact 
functional form of the decrease because of the subsequent 
increase in k and the unusual peak at —5.5 K. Neverthe- 
less, if we confine ourselves to T < 5 K, the data can 
be analyzed; assuming that k s diminishes below Tsp, 
we can assume n is only due to phonons for T <C Tsp. 
Using our heat capacity result, we can divide n by the 
lattice specific heat, c p = /3T 3 . The resulting diffusivity 
in units of cm 2 /sec is displayed in the inset to Fig.l, plot- 
ted against T -1 . Remember that we expect the plotted 
D to represent D p only below —5 K (T _1 >0.2 K _1 ). In 
this temperature range it is clear that D — T _1 , a tem- 
perature dependence commonly attributed to a phonon 
scattering rate due to planar defects such as the strain 
field of dislocations |p2| . In fact, transmission electron 
microscope investigation of our crystal found planar de- 
fects whose average distance is about 10 /im p3]] . 

Based on the low-temperatures result discussed above, 
we can assume a model in which n is dominated by 



phonons below — 12 K and exhibiting a competition be- 
tween two scattering mechanisms that cause the peak 
at —5.5 K. If we assume phonons are scattered both 
by defects and by spin excitations, we can write D p = 
[AT + £(A/T) r e- A / T ]-\ where A and B are constants 
and r is a power smaller than n (as introduced above). 
Here, A can be determined from the low-temperature be- 
havior (inset of Fig.l). Dividing n by c p = /3T 3 again and 
assuming 1 < r < 2 (r has to be in that range to accom- 
modate the spin-excitations density of states as found 
earlier and to give a maximum to the formula), we esti- 
mate B — 10 _3 sec/cm 2 . It is easy to calculate now the 
phonon mean free path due to phonon-spinon interac- 
tion, £ p - s . Using the estimates for B and r and a sound 
velocity of — 5 x 10 5 cm/sec, we find l p - s (T = T S p) —10 
/im; this value is compatible with the phonon mean free 
path for planar-defect scattering, £ p -d, estimated from 
the inset of Fig.l (£ p -d~20 /im at 5 K), which agrees 
with the actual planar-defect distance of —10 /im. 

Figure 2 shows k(T) in the field range - 16 T. One 
striking feature here is the strong suppression of the 5.5-K 
peak with magnetic field; k changes by more than a factor 
of 3 at 5.5 K, a big magneto-thermal effect. This strong 
effect of the field on the peak is another confirmation that 
its origin is related to the magnetic excitations in the SP 
phase. Note that the position of the peak does not shift 
much as a function of temperature for magnetic fields be- 
low 12 T, similar to the gap that does not change much. 
The peak becomes a shoulder at 14 T and a new low tem- 
peratures upturn below 4 K emerges at higher fields. At 
the same time, the higher temperatures data is almost 
unchanged with the broad maximum at —22 K appear- 
ing for all fields. The transition itself is clearly shifted 
towards lower temperatures exhibiting similar change in 
slope for all fields below —12 T. In the inset to Fig. 2 we 
plot Kj as a function of magnetic field at 4.2 K to demon- 
strate the behavior of n near the peak. Besides the sharp 
drop in n at —12.5 T, a hysteresis is found upon cycling 
the magnetic field. This effect is a clear manifestation of 
the first-order transition to the I phase at high magnetic 
field. 

Lorenz et at. recently reported specific heat measure- 
ments of CuGeOs in 16 T field (l7| . The low-temperature 
(T < 6 K) behavior of their 16 T specific heat data is 
consistent with c = ((3 ph + Pm ag )T 3 with (3 mag -1.4 
mJ/mole-K 4 . Since in this high field and low temper- 
ature both phonon and magnetic contributions give a T 3 
dependence to c, one can divide k by T 3 to search for dif- 
ferent temperature dependences of the scattering times. 
Again, examining k(T) in 14 and 16 T, we find that the 
diffusivity is proportional to 1/T below T —4 K. Since 
k at low temperatures is greatly suppressed and the dif- 
fusivity has still the 1/T dependence, we may conclude 
that phonons are strongly scattered by solitons in the I 
phase. Also, the smallness of k in the I phase compared 
to the SP phase suggest that solitons do not carry much 
heat. Note that the solitons are simultaneously magnetic 
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and structural excitations: they carry spin 1/2 and are 
domain walls in the dimerized lattice. 

Finally, we briefly discuss the increase in k in the I 
phase from 14 to 16 T. The T 3 dependence of the mag- 
netic specific heat in the I phase was explained by Bhat- 
tacharjee et al. [|I| as the specific heat of the gapless 
phason modes. The increase of this phason with increas- 
ing field can be the cause of the increase in ft, although 
the magnetic field dependence of the phason has not been 
calculated. Another possibility is that the solitons them- 
selves carry heat and the soliton contribution is propor- 
tional to the soliton population. However, since there is 
no theory for either soliton transport or phonon-soliton 
scattering, we cannot check the consistency of the soliton 
scenario. 

To summarize, our thermal conductivity measurement 
of CuGeOs revealed four characteristic features: (a) a 
broad peak centered at ~22 K above the SP transition, 
(b) downward kink at Tsp as temperature is lowered, (c) 
unusual peak centered at ~5.5 K below T^p, (d) strong 
suppression of the 5.5-K peak with high magnetic fields. 
These observations are consistently understood by con- 
sidering both phonon and spinon heat transports, which 
are analyzed in combination with specific heat data. In 
particular, the unusual 5.5-K peak is a consequence of 
phonons scattered both by defects and by spin excita- 
tions. 
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FIG. 1. Thermal conductivity and specific heat of 
CuGeOs single crystal in H=0. Inset: Low temperatures 
diffusivity extracted by dividing k by the fitted lattice 
specific heat. 

FIG. 2. A set of k(T) in magnetic fields up to 16 T. 
Inset: Kj as a function of magnetic field at 4.2 K. Note 
the hysteresis at the SP to I transition at ~12.5 T. 
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